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In  conjunction  with  adjoint  variational 
principles,  two  classical  problems  of  elastic 
stability  under  nonconservative  forces  and 
under  the  influence  of  internal  and  external 
damping  are  studied  using  the  finite  element 
technique.  The  solutions  are  formulated  in 
the  framework  of  Rayleigh-Ritz  method.  It 
is  demonstrated  that  this  approach  is  very 
effective  for  solutions  of  complicated  non¬ 
conservative  stability  problems.  The 
destabilizing  effect  of  Ziegler  due  to 
internal  damping,  the  extreme  sensitive 
nature  of  stabilizing  or  destabilizing 
effects  for  very  small  internal  or  external 
damping  parameters  and  for  other  ranges  of 
such  effects  of  practical  interest  have 
been  obtained  for  these  two  problems  and 
are  presented  in  several  sets  of  curves. 


Columns 

Stability 

Vibrations 

Damping 

Finite  Element 

Nonadjoint  Variational 
Principles 


1 


TABLE  OF  CONTENTS 


Page 

ABSTRACT  ! 

INTRODUCTION  3 

STATEMENT  OF  THE  PROBLEM  AND  THE  BASIS  OF  SOLUTIONS  5 

FINITE  ELEMENT  FORMULATIONS  9 

METHODS  OF  COMPUTATION  14 

RESULTS  AND  DISCUSSION  16 

APPENDIX  .  18 

REFERENCES  21 

DD  Form  1473  (Document  Control  Data  -  R§D) 

LIST  OF  CAPTIONS 

Fig.  1  Problem  Configurations  and  the  Finite  Element  Idealization  25 

Fig.  2  Frequency  vs.  Load  for  the  Two  Lowest  Modes  of  Leipholz' 

and  Beck's  Problems  Without  Damping  26 

Fig.  3  Leipholz'  Problem:  Qcr/1f2  vs.  o  27 

Fig.  4  Leipholz'  Problem:  Q^/w2  vs.  Small  n  (n^0.5)  28 

Fig.  5  Leipholz'  Problem:  Q^/*2  vs.  B  29 

Fig.  6  Leipholz;  Problem:  Q^/"2  vs.  Small  B  (B^0.5)  30 

Fig.  7  Beck's  Problem:  Q^/tr2  vs.  n  31 

Fig.  8  Beck's  Problem:  Q^/*2  vs.  Small  n  (n^0.02)  32 

Fig.  9  Beck's  Problem:  Q^/ff2  vs.  6  33 

Fig.  10  Beck's  Problem:  Q^/*2  vs.  Small  6  (Bj«5.0)  34 


2 


1.  INTRODUCTION 


The  problems  of  elastic  stability  under  nonconservative  forces  have 
attracted  many  researchers  in  the  last  two  decades.  A  review  paper  by 
G.  Herrmann  [1]  contains  a  comprehensive  bibliography  of  work  performed 
before  1967.  The  practical  significance  of  this  class  of  problems  has 
become  increasingly  noted  with  the  advancement  of  modern  technology.  A 
flexible  missile  under  thrust  and  a  structural  p.-rt  of  an  aircraft  under 
aerodynamic  forces  are  obvious  examples.  Approximate  methods  are  usually 
required  for  the  solutions  of  such  problems.  Galerkin's  method  was  first 
used  by  H.  Leipholz  [2]  to  obtain  solutions  for  a  nonconservative  system. 
Because  boundary  value  problems  in  the  theory  of  nonconservative  elastic 
stability  are  nonself-adjoint,  by  virtue  of  the  fact  that  nonconservative 
forces  do  not  possess  potentials,  no  complete  functional  exists  for  the 
classical  form  of  Hamilton's  principle.  M.  Levinson  [7]  points  out  later 
that  the  conventional  Hamilton's  principle  with  nonconservative  forces 
can  be  written  as  a  well-posed  variational  principle  with  some  constraint 
conditions  and  the  Ritz  method  can  thus  be  employed,  formally  at  least, 
for  the  solution  of  these  problems.  A  more  recent  trend  has  been  the 
interest  in  the  establishment  of  well-posed  variational  principles  with¬ 
out  any  constraint  conditions  by  the  introduction  of  the  adjoint  systems 
[4-10]. 

Based  on  the  variational  principle  suggested  by  Levinson,  the  finite 
element  method  has  been  employed  for  solutions  of  several  nonconservative 
stability  problems  without  damping  [11,12,13]-  In  contrast,  the  present 
work  is  an  application  of  the  variational  principle  of  Leipholz'  type  [5]. 
The  effects  due  to  both  the  internal  and  external  damping  are  also  included 
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in  this  study.  The  finite  element  method  has  been  shown  to  be  of  the  Ritz 
type  of  approximations  in  the  linear  theory  of  elasticity  [14,15,16].  A 
variational  principle  is  therefore  a  natural  basis  for  its  formulations. 

The  effect  of  damping  on  various  nonconservative  stability  problems  has 
been  studied  by  many  investigators  [9,10,17-34].  The  difficulty  arises 
not  only  from  the  nonself-adjoint  nature  of  the  differential  equations 
but  also  from  the  fact  that  the  eigenvalues  involved  are  complex  numbers 
in  general. 

Two  classical  problems  with  both  internal  and  external  damping  terms 
are  considered  in  this  work,  i.e.,  the  Leipholz'  problem  [2]  (a  cantilevered 
column  with  uniformly  distributed  force  tangential  along  the  column's  length) 
and  the  Beck's  problem  [35]  (a  cantilevered  column  subject  to  a  follower 
force  at  its  free  end) .  The  internal  damping  is  assumed  to  result  from  a 
Kelvin-Voigt  type  of  material  and  the  external  damping,  due  to  a  dissipative 
force  proportional  to  the  velocity  at  a  point  in  the  column.  In  the  next 
section,  the  differential  equations  and  boundary  conditions  of  the  problems 
and  their  ad joints  are  defined.  The  derivation  of  the  respective  adjoint 
variational  principles  is  also  outlined  here.  In  Section  3,  we  give  the 
finite  element  formulations  that  lead  to  the  matrix  eigen-equations.  Final 
results  and  discussion  on  the  numerical  methods  employed  in  this  work  are 
presented  in  Sections  4  and  5. 
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2*  -SJ.ATEMENT  OF  THE  PROBLEM  AND  THE  BASIS  OF  SOLUTTflNfi 

The  configurations  and  loads  of  Leipholz'  problem  and  of  Beck's  are  shown 
in  Figures  1(a)  and  1(b).  In  the  absence  of  damping,  the  differential  equa- 
tion  for  Leipholz  problem  [2]  is 


El 


.4- 

d  u 


.2- 

9  u 


+  q(*-x)  *  pa  iJL  = 

3x 


ax  3x^ 

and  for  Beck’s  problem  [33], 


af2 


o. 


*4-  .2— 

PT  9  U  _  3  n 

EI  "T3T  +  P  +  PA 
3x  3x2  3t 


=  0, 


(2) 


with  boundary  conditions  for  both  cases: 


x=0 


3u 

ax 


=  0 


and 


a2u 

3x2 


x=0 


3_ 
a  u 

ax3 


=  o 


(3) 


x=t 


where  u  -  uix.t)  is  the  deflection  of  the  column  from  its  undisturbed  postion; 
x.t  are  the  spatial  and  time  variables,  respectively;  E,  the  Young's  modulus, 
o.  the  density;  1,  the  second  moment  and  A,  the  area  of  the  cross-section.  The 
letter  q  denotes  the  force  per  unit  length  in  the  Leipholz'  problem  and  P  the 
concentrated  follower  force  in  Beck’s  problem. 

To  include  damping  effects,  a  dissipative  for  F4,  such  as  the  air  resistance, 
is  assumed  to  be  proportional  t-  the  velocity  of  the  column  at  a  point,  i.e., 


Fd  ’  -V3 

1  3t 


(4) 
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where  8j  is  the  external  imping  constant.  The  material  is  assumed  to  be 
of  a  linear  viscoelastic  solid  of  a  Kelvin-/oigt  type.  The  one  dimensional 
stress-strain  relation  can  be  written  as 

de 

a  =  E  e  +  E*  —  (5) 

dt 


where  o,e  are  the  one-dimensional  stress  and  strain  respectively;  and  E*, 
the  viscosity  of  the  material. 


The  equations  of  motion,  including  damping  terms,  are  then  [9] 


34u 

El  — _ tE*I 

37* 


*5— 
3  u 


3)c4dt 


3Zu  3u 

♦  q(£-x)  —  ♦  B,  —  +  PA 
3X2  at 


0 


for  Leipholz'  problem  and 

34u  3^u  3^u  3u  3^u 

El  +  E*I  +  P  —  +  8,  —  ♦  PA  =  0 

3x  3x  3t  3xZ  3t  3t 


(6) 


(7) 


for  Beck's  problem. 

It  will  be  convenient  to  use  nondimensional  quantities  for  the  analysis. 


Thus, 

let 
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El 

where 

the 

constant 

c  has  the  dimension  of 

time  and  all 

other  quantities  in 

Eqs.  (8)  are  dimensionless.  The  nondimensional ized  load  parameters  are 


Q  =  - 

El 
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and 
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for  the  adjoint  to  the  Leipholz  problem,  and 


r™  -  rtvIV  ♦  Qv"  +  v  -  8  v  =  0 


(15) 


V  |  =  v'  !  =0 

x=0  x=0 

[v''  -  n  v' *  ♦  Qv]  =  [v1 "  -  nv' ' '  +  Qv*]  -  0 

x=l  x=l 


(16) 


for  the  adjoint  to  Beck's  problem,  where  v  denotes  the  adjoint  field  variable. 

For  Leipholz5  problem  and  its  adjoint,  one  can  consider  Eqs.  (10)  and  (13) 
as  the  dynamic  equations  of  equilibrium  of  the  original  and  the  adjoint  system. 
Introducing  virtual  displacements  <Hi  and  $v  which  satisfy  the  respective  end 
conditions,  it  is  then  possible  to  generalize  the  principle-  o'.  .’vrtual  work 
such  that 

[uIV  +  n  uIV  ♦  Q(l-x)u"  ♦  8u  +  u]6v 


♦  [vIV  -  n  vIV  ♦  {Q(l-x)v)"  ♦  v  -  Bv]$u  =  0. 


(17) 


Integrating  Eq.  (17)  over  the  interval  (0,1)  and  between  the  two  time  limits  tj 
and  t2>  and  imposing  end  conditions  of  Eqs.  (11)  and  (14),  a  generalized  varia¬ 
tional  principle  can  be  written  [9]  as  the  following 


6  /  2  L  dt  *  0, 

.  ,  h 

with 

L  =  T  -  V  -  Vd, 
and 

T  =  i  /  u  v  dx, 

2  o 

1 

V  =  1  /  (u"v"  -  Q[(l-x)v]'u’)dx, 
2  o 


(18) 

(19) 

(20) 
(21) 
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(22) 


1 

V .  =  ~  /  ( n  (v  ’ '  u '  *  -  u"v")  +  6 (vu  -  uv)}dx, 

4  o 

can  be  referred  to  as  the  generalized  kinetic  energy,  potential  energy  and 
dissipative  potential  resepctively. 

Similar  results  can  be  obtained  for  Beck's  problem  with  the  expression 
of  V  'in  Eq.  (21)  replaced  by  that  of  Eq.  (23).  Thus 

1  1  i 

V  =  i  /  [u"v"  -  Q  u'v’Jdx  +  i  Q  u' (l,t)v(l,t)  (23) 

2  o  2 

for  Beck's  problem  and  its  adjoint. 

The  finite  element  formulations  are  based  on  the  variational  equation  of 
(18),  where  L  is  clearly  a  complete  functional.  In  comparison,  the  work  by 
Barsoum  [11]  and  that  by  Mote  [12,13]  are  based  on  an  extended  Hamilton's 
principle  of  the  form 


b 


<5  /  2  L,  dt  +  /  2  6W  dt  =  0  (24) 

tl  tl 

where  is  a  complete  functional  but  6Wlt  which  is  the  virtual  work  done  by 
nonconservative  forces,  cannot  be  express °d  as  the  variation  of  a  path- 
independent  functional.  Using  a  different  approach,  Anderson  [36]  has  shovn 
that  the  adjoint  variational  principle  (18)  can  lead  to  more  accurate 
numerical  solutions  than  does  the  variational  expression  of  (24) . 

3.  FINITE  ELEMENT  FORMULATIONS 

In  the  application  of  the  finite  element  method,  the  column  being  analyzed 
is  divided  into  L  segments  (elements)  as  shown  in  Figure  1(c).  In  conjunction 
with  the  variational  principle  of  Eq.  (18),  the  integrations  in  Eqs.  (20),  (21) 
and  (22)  are  replaced  by  summations  of  the  quantities  in  all  the  elements.  (The 
following  derivation  is  carried  out  for  Leipholz'  problem.  Similar  outline  of 
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the  procedure  for  Beck's  problem  will  not  be  repeated  here  [37].)  Thus, 


T  =  V  if1  a(i5v(i)dx 


Jl  *  l 


i-1 


V  = 


j  j  j1  ju"^"^  -  Q[(l-x)v(l)]'  u'^Jdx 
i-i 


(25) 


(26) 


and 


Vd  = 


L  £ 

^  -  r  1  [n(v",fl^"Cl^-  u',(l)v"(l))  ♦  B(vu  -  uv)]dx  (27) 

i-i  2  ij.j 


where  l  =  0,  £L  *  1;  and  the  superscript  (i)  denotes  the  quantities  in  the 
i-th  element. 

We  effect  the  following  change  of  independent  variable.  Let 


C  =  — -  (x  -  ^.j)  *  L(x 


i-1 


) 


(28) 


1  1-1 


In  Eq.  (28)  and  in  the  sequel,  we  have  assumed  that  all  the  elements  are  of  the 
same  length.  Consequently, 

d? 


dx 


=  L, 


d2C 


37 


=  0 


(29) 


For  some  function  f(x), 

df  ,df  d2f 


=  L 

dx  d?  dx 

where  the  notation 

f(x)  -  f[x(?)]  =  f(C) 

has  been  adapted  for  simplicity. 


d2f 

d?2 


(30) 
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Using  the  new  independent  variable  5,  Eqs.  (25),  (26)  and  (27)  become  the 
following: 


1  y  r*  •(i)*(i),r 

T  =  —  1  )  u  v  d? 

2L  i=l  o 

V  =  -  l  /  (L3u" 

2  i=l  o 

+  Q[u,(iV1J-  (L-i+l-e)u,(iV(l)])dt 


and 


Vd  =  ~  l  J  [nL3(v"(i)u*,fi)-  U"fi)v"(i)) 


i=l  o 


(25a) 


(26a) 


+  £  (v(l)u^-  u(i)v(i))]d?.  (27a) 

L 

We  shall  use  a  polynomial  to  approximate  the  displacement  field  in  an  element. 
The  displacement  and  its  first  derivative  (slope)  will  be  required  to  be  con¬ 
tinuous  between  the  two  adjacent  elements.  Hence  there  are  four  degrees  of 
freedom  for  each  element  and  the  polynomial  shall  be  cubic.  In  terms  of  a 
Rayleigh-Ritz  approximation,  the  coordinate  function  is  chosen  to  be  a  piece- 
wise  analytic  function. 

It  is  convenient  to  let 

u(l)(?,t)  =  aT(?)  X(t)  =  aT(i;)  U(l)eH 
and 

v(i)(C,t)  =  aT(C)  Y(t)  *  aT(?)  V^e** 
where  X  is  the  usual  eigenvalue  parameter, 

aT(0  =  {ax(c)  a2(c)  a3U)  a4(c)> 
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is  the  displacement  function  vector  and 


UT(iJ  =  u2(1)  usC» 

yT(i)  _  |y^(i)  v  (i)  y^(i) 

are  the  generalized  displacement  vectors  in  the  i-th  elementi  and  the  superscript 

T  denotes  the  transpose  of  a  matrix  (or  a  vector).  When  we /take  U2^  to 

» 

be  the  displacement  and  the  slope  respectively  at  the  left  jsnd  (see  Figure  1(c)) 
of  i-th  element  and  take  Uj^,  U4^  to  be  the  same  at  the  right  end  (same  for 
the  vector  V^),  it  is  easy  to  see  that  a(S)  must  take  the  following  form: 

T 

a  (?)  =  (aj(?)  a2(?)  a3(C)  a4(C)>  / 

i 

2  X  9  r  7  X:  7  x  C34) 

=  (1-3?  +  2?  C-2<,  +5  35  -25  •  -5  +53). 

Using  Eqs.  (28) -(34)  in  the  generalized  energy  expressions  Eqs.  (25a),  (26a)  ana 

(27a),  and  carrying  out  the  variation  of  Eq.  (18),  we  obtain  the  following 

variational  equation  for  the  descrete  system: 


v(i» 

4 


(33) 


t., 

6  /  2  L  dt 

=  -J2  l  e2At{«UTCiJ{X2-  -  A[nL3C  ♦  8  f 

2  t  i*l  L  1  -  L 

1 

♦  L3C  ♦  Q[H  ♦  £  -  (L  -  i+l)Bj)V(l) 


♦  6VT(1){X2  -  ♦  A[nL3C  ♦  8  -  ] 

L  L 

♦  L3C  Q[E  +  FT  -  (L  -  i-*l)BjU(i)}}dt  *  0  (35) 

where  A,  IB,  £  and  £  are  constant  matrices  whose  definition  and  numerical  values 
are  given  in  the  Appendix.  It  is  noted  in  the  above  mentioned  matrices  only  £ 
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is  not  symmetric.  In  Eq.  (35),  6U^  and  6V^  are  independent  of  each 

,(i) 


other.  Since  we  are  interested  only  in  the  solution  of  IP1^,  we  can  set 

.i  =  1.....L.  (36) 


<5UCl5  =  0, 


Eq.  (35)  now  becomes 
L 


tj  i=l 


fZ2  r  2^t  -,,T(i),.2  —  r  3_ 

)  l  e  o/  (X  —+X[nL£+6  —  ] 

L 


+  L3C  +  Q[E  +  FT  -  (L-i+l)iy  HJ^dt 


=  0 


(37) 


Up  to  this  point,  the  matrices  involved  are  referred  to  the  individual  elements 
with  a  size  of  4x4.  Next  we  must  form  the  overall  matrix  equation  and  eliminate 
the  redundant  unknowns.  To  accomplish  this,  the  boundary  conditions  and  the 
continuity  conditions  must  be  applied.  In  Leipholz'  problem,  the  boundary 
conditions  for  the  original  and  the  adjoint  problem  has  the  same  form  as  the 
following  (see  Reference  [37]  for  a  discussion  on  the  adjoint  boundary  conditions 
in  Beck's  problem.) 


r  (i)i 
V 

i - 

© 

i  ■■■■-. 

< 

► 

* 

UC1). 

L  °. 

Vu' 

i 

1 

'  u/L> 

► 

s 

- 

U™ 

0 

1 

.  U2a’. 

(38) 


The  continuity  conditions  are  also  the  same  for  the  original  and  the  adjoint 
problems: 


u3<1-1)  =  UjW 


u  Ci-D 

4 


u2(i) 


2,3, .... ,L. 


(39) 
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Introducing  new  vectors  with  independent  elements 


UT  =  (U,^'  U il  U  <2> . 

3  4  3  4  3  4  (40) 

VT  =  {V  v  v,(2)  . v  (L-l)  u  (L-i) } 

3  4  3  4  3  4 

and  using  the  conditions  of  Eqs.  (38)  and  (39),  Eq.  (37)  can  be  rewritten  as 


»2At 


KUdt  =  0 


(40) 


where 

K  =  A2  R  +  AS  +  T  (41) 

is  a  constant  matrix  of  2 (L-l )  x  2(L-1)  and  R,  S  and  T  are  also  given  in  the 
Appendix. 

Since  now  the  vector  6V  in  Eq.  (40)  is  arbitrary  and  all  its  elements  are 
independent,  we  obtain  the  final  matrix  eigen-equation 


K  U 

*  (A2  R  +  AS  +  T)U  =  0  (42) 

which  is  to  be  solved  next. 


4.  METHODS  OF  COMPUTATION 

Because  dissipative  systems  possess  complex  eigenvalues  in  general,  let 
A  =  Aj  +  iAj,  i  =  /T  (43) 

denote  the  eigenvalue  of  the  system,  where  A^  and  A2  ?re  real  numbers.  When 
Aj  is  positive  and  A2  t  0,  the  system  loses  stability  by  flutter,  i.e.  oscillations 
with  increasing  amplitude.  For  each  set  of  damping  parameters  8  and  n,  the 
critical  load  Qcr  is  obtained  when  Aj  has  just  changed  sign  from  negative  to 
positive.  In  this  study,  two  methods  were  used  in  extracting  the  eigenvalues. 
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They  are  described  briefly  here. 


A.  Standard  QR  Algorithm  [38,39]  - 

Three  steps  are  required  in  using  this  approach: 

1.  To  transform  Eq.  (55)  into  the  standard  eigenvalue  equation 

(M  -  xI)U  =  0  (44) 

where 


r  *  i 


and  0,  I  are  zero  and  unit  matrices,  respectively. 

2.  To  transform  the  matrix  M  into  Hessenberg's  form. 

3.  To  extract  all  the  eigenvalues  from  a  Hessenberg's  matrix  using  the 
QR  iteration.  Steps  (2)  and  (3)  are  accomplished  by  two  Standard  Subroutines 
provided  in  the  IBM  System/360  Scientific  Subroutine  Package,  i.e,  HSBG  and 
ATEIG,  respectively. 

B.  Rosenbrock's  Method  of  Iteration  [40]  - 

This  is  a  trial  and  error  method  using  the  route  along  which  the  function 
decreases  most  rapidly. 

In  terms  of  computer  time.  Method  B  is  less  efficient  than  Method  A.  In 
applying  Method  B,  there  is  a  great  tendency  to  converge  to  the  wrong  eigen¬ 
value.  However,  once  located  correctly,  the  eigenvalue  can  be  determined  with 
more  accuracy  than  Method  A. 
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S.  RESULTS  AND  DISCUSSION 


Calculation  for  critical  loads  Qcr  are  performed  for  the  range  of  damping 
parameters  which  is  both  of  practical  interest  and  sensitive  in  stabilizing 
or  destabilizing  effects.  Results  are  shown  in  Figures  2  through  10. 

In  Figure  2,  the  two  lowest  branches  of  the  frequency  curves  are  shown 
for  both  the  Leipholz'  and  the  Beck's  problem  without  either  internal  or 
external  damping.  For  small  load  Q,  A  is  a  purely  imaginary  number  (Aj  =  0, 

A2  =  in  the  notation  used  in  the  previous  section.)  As  the  value  of  Q 
increases,  eventually  the  two  lowest  branches  coincide  and  is  reached, 
beyond  which  Aj  becomes  positive  and  the  column  will  fail  by  flutter.  In 
Figure  3  through  Figure  10  are  presented  Q/**  vs.  the  damping  parameter  B  or 
n.  These  curves  have  shown  that  for  the  range  0<B<2.0  and  0<n<0.2,  the  damping 
effect  is  extremely  sensitive.  In  this  range,  the  internal  damping  has  a 
destabilizing  effect  as  can  be  seen  from  Figures  3,  4,  7  and  8.  The  internal 
damping  can  also  have  stabilizing  effect  in  certain  ranges,  e.g.  the  range  of 
0<6<5.0  and  0.1<n<0.5  for  Leipholz'  problem  as  shown  in  Figure  4.  On  the  other 
hand,  the  external  damping  always  has  a  stabilizing  effect  which  is  very  sensi¬ 
tive  when  both  parameters  are  small.  This  is  seen  in  Figures  5,  6,  9  and  10. 

It  may  be  of  some  interest  to  note  in  Figure  4  for  Leipholz  problem  that 
when  6=0  and  n  approaches  zero,  the  limit  of  Qgj.  is  2.32*  while  in  Figure  5, 
when  <i  has  been  set  to  zero  to  begin  with,  Q  has  a  value  of  3.93*2.  For 

cr 

Beck’s  problem  similar  situation  is  observed  in  Figures  8  and  10  with  Qcr|^~o  * 
l.lOn^  an(j  B=0  _  2.03*^.  These  are  the  cases  when  the  solution  of  the  limit 
of  the  problem  is  not  the  same  as  the  limit  of  the  solution  of  the  original 
problem  —  a  phenomena  known  as  the  destabilizing  effect  of  damping  first  noted 
by  Ziegler  [41]. 
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A  few  words  about  the  convergence  of  the  data  are  in  order  here.  Both 
Barsoum  [11]  and  Mote  [12]  have  shown  that  for  similar  problems  and  using 
the  same  displacement  functions,  the  finite  element  solutions  using  eight 
segements  converge  to  within  1%  of  the  known  exact  solutions.  Eight  is  the 
number  used  in  this  analysis.  The  solutions  of  Beck's  problem  presented 
here  for  external  damping  alone  show  excellent  agreement  with  the  data 
obtained  by  Plaut  and  Infante  [27] .  In  comparing  with  the  solutions  for 
Leipholz'  problem  without  damping  obtained  by  McGill  [42]  using  a  two-term 
Galerkin's  procedure  and  those  obtained  by  Anderson  and  Walter  for  both 
Leipholz'  and  Beck's  problem  with  both  using  a  three-term  Ritz  approxima¬ 
tion,  the  data  here  show  general  agreement.  Good  agreement  has  also  been 
established  between  the  data  obtained  here  and  those  by  Dzydlo  and  Solarz 
[28]  from  the  numerical  calculations  based  on  the  exact  frequency  equation. 
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appendix 


Definitions  and  values  of  some  matrices  appeared  in  Section  3: 

r 1  T 

I  a(C)  a  (Qd5 


1  4 

5  +  7  (SYMMETRIC) 


1  (A-l) 

/  a’(c)  aT(c)d  t, 


(SYMMETRIC) 


-6 


1  -  3  ) 


l  i 
T  "  I 


•!6 


/  a"(?)  aT'»(?)dc 

o 


I  5  a'(c)  aT'(?)d? 
o 


(SYMMETRIC) 

12 

-6  4 


(!-!♦  1 

4  S  W 


-  1 


5 

0 


1  .  1 
5  6 

-  _L 

10 

_  1 
6U 


(SYMMETRIC) 


3 

S 

0 
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(A-2) 
(A- 3) 


(A-4) 


/  a*  (?)  aT(?)d? 
o 


R  = 

rM: 


l 

_  1 

1 

10 

2 

10 

I  +  i 

1 

1 

2  6 

"  10 

60 

1 

1 

1 

10 

2 

’  10 

_1_ 

_1_ 

o 

60 

10 

1  rC2) 

R13 

(2) 

R14 

0 

1  r(2) 

K23 

r(2) 

R24 

0 

r (2)  r(2)+r(3)  r(2)+r(3) 

32  33  11  34  K12 

42  43  21  44  22 


;a-2)  (L-1)  (L-2)  (L-l)  (L- 

R33  R11  R34  +R12  R13 

RfL-2)+R(L-l)  (L-2)  (L-l)  (L- 

R43  +R21  R44  +r22  R23 


r(L-1) 

k31 


,(L-1) 


(a-l)  d(L-1) -d 


R$"J)  rJ^+R 


R<»  .  *g>  *  RC«  *  RC« 


RU  .  R#  *  R{«  .  RgJ  .  Rt«  .  R&> 

*fi»  *  *$>  *  *  «g>  *  «ft>  *  «g> 

r|«  .  r||>  *  r£>  *  Rg)  ♦  Rft)  *  R«  *  R<J1  ♦  R$>  *  rJJ) 


The  expressions  for  S  and  T  are  obtained  by  replacing  the  letter  R  with  the 
letter  S  and  T  respectively  in  Eqs.  (A-6)  and  (A-7)  wherever  R  appears. 
Furthermore,  R(l),  S(i>  and  are  defined  below. 


rW  = 

L 


(A- 8) 


SW  =  n  L3  C  +  6  - 

L 


(A-9) 


=  L3  C  +  Q[E  +  rT  -  (L  -  i  + 


(A-10) 


where  i  *  1,2, .... ,L. 
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(c)  FINITE  ELEMENT  MODEL 


rc  I. 


PROBLEM  CONFIGURATIONS  AND  THE 
FINITE  ELEMENT  IDEALIZATION 


FREQUENCY  VS.  LOAD  FOR  THE  TWO  LOWEST  MODES  OF 
LEIPHOLZ'  AND  BECK'S  PROBLEM  WITHOUT  DAMPING 


Qcr/it2 


igure  5. 
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LEIPHOLZ1  PROBLEM:  QcR/r2  VS.  0 


Figure  7 

BECK'S  PROBLEM1  Qcr/it2  VS.  i) 
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